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Abstract. To date, classical mass renormalization has been invoked in all tube-based 
self-force calculations, thus following the method introduced in Dirac's 1938 calculation 
of the electromagnetic self-force for the classical radiating electron. In this paper a new 
tube method is described that does not rely on a mass renormalization procedure. As 
a result, exact self- force calculations become possible for classical radiating systems of 
finite size. A new derivation of the Lorentz-Dirac equation is given and the relationship 
between the new tube method and the classical mass renormalization procedure is 
explained. It is expected that a similar tube method could be used to obtain rigorous 
results in the gravitational self- force problem. 
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1. Introduction 

The classical mass renormalization procedure was introduced 70 years ago by Dirac pQ 
and has since featured in all tube-based calculations of self-forces on radiating particles. 
This includes not only calculations involving the electromagnetic self-force on a radiating 
charge in Minkowski spacetime, as exemplified by [T] [T7] , but also electromagnetic self- 
force calculations for a charged particle in a curved background spacetime [TBI [19] . 
and the tube-based approach to the gravitational self-force problem, as taken in [20j 
and reviewed in [21]— [24J . Nevertheless, the mathematical status of the classical 
mass renormalization procedure has remained unclear and justification for its use has 
principally been that it has given results that agree with other self-force methods. 
As recently emphasised in [25], this is particularly true of the gravitational self-force 
problem for which the concept of a point particle makes even less sense than in the 
electromagnetic case. 
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This paper describes a tube method for determining the self-force on a radiating 
particle or system of finite non-zero size. Dirac's tube method is extended so that 
the calculations go through with a worldtube that has, at all stages of the calculation, 
a non-zero radius. The physical picture that emerges is very different to that which 
accompanies the classical mass renormalization procedure. In retrospect, one finds that 
the renormalization procedure is most easily understood as a short-cut that involves 
assuming, rather than deriving, a certain intermediate result (our equation (22)) that 
is required for a proper derivation of the self-force. 

The new tube method is presented here for the problem of electromagnetic self- 
force in Minkowski spacetime and, as an example, a derivation of the Lorentz-Dirac 
equation is given that makes no mention of mass renormalization. Since the method 
is based on quite general considerations of energy-momentum conservation, the same 
approach ought to be applicable in any radiation dynamics problem for which a well 
defined self-force exists. In particular, it seems likely that a similar tube method could 
be used to obtain rigorous results in the gravitational self-force problem. 



1.1. Background 

Dirac's tube method for calculating the self-force on a radiating particle was used in 
PQ to derive what is now known as the Lorentz-Dirac equation of motion for a classical 
radiating particle of charge q and mass m in a background electromagnetic field F, 

mcz a = qF* p iP + I (8<% + z a zp) z* , (1) 

where the worldline has equation x = z(s) in Minkowski coordinates, overdots denote 
differentiation with respect to s, and the parameterization is such that z a z a = — 1. 

The calculation that led Dirac to (II]) can be described as follows. One starts 
by deriving a balance relation that follows from energy-momentum conservation, as 
applied to the electromagnetic field within a worldtube containing the radiating particle. 
The balance relation expresses the fact that the rate of change of 4-momentum flowing 
within the worldtube must equal the rate at which electromagnetic energy- momentum 
is radiated through the wall of the tube. The derivation of Dirac's balance relation is 
lengthy but straightforward. It is easily found (for a given worldtube) using computer 
algebra and takes the form of a power series in the tube radius r, which is taken to be 
small compared with the radius of curvature of the particle worldline (the inverse of the 
particle acceleration). 

The worldtube exists only to facilitate the calculations, so the tube radius can not 
enter into the sought-after equation of motion. Dirac's approach for eliminating the 
tube radius was to let r — > 0, in which case only quantities determined by the particle 
worldline can appear in the final result. The difficulty with this idea is that because the 
particle must remain enclosed within the worldtube, it must then be modelled as being 
point-like. The electromagnetic field is then singular on the particle worldline, causing 
the required r — > limit to be neither defined nor physically sensible. In fact, to obtain 
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the desired result ([T]) using this approach, one must formally subtract out the divergent 
Coulomb field energy of the point charge by assigning the particle an infinite negative 
bare mass — the procedure that has come to be known as classical mass renormalization. 

To avoid the absurdities of classical mass renormalization, one needs a way to 
eliminate the tube radius from Dirac's balance relation without taking r — > 0. Moreover, 
simply taking r to be small and neglecting O(r) terms does not quite work, as this would 
result in a finite mass renormalization that is r-dependent. In effect, one would be saying 
that the field-energy contribution to the mass of the particle somehow depended on the 
radius of an imaginary tube that was introduced purely as a mathematical convenience. 

In this paper we shall see how to cancel all of the r-dependent terms from the 
Dirac balance relation for any classical radiating system. This is done by making better 
use of the information available in the problem, including the fact that the background 
electromagnetic field must satisfy the source-free Maxwell equations throughout the 
spacetime region occupied by the radiating particle or system. The worldtube used in 
the calculations will have a non-zero radius and there will be no need to assume any 
mass renormalization, finite or otherwise. 



1.2. Organisation of the paper 

The main ideas and results are presented in Section [2j where the tube method is 
described with reference to the Lorentz-Dirac example, that is, for the problem of 
determining the self-force for a charged particle in Minkowski spacetime. As far as 
possible, all distracting calculations have been deferred to later sections. In particular, 
the geometry of the worldtube is dealt with in Section [3] and the electromagnetic field 
calculations appear in Section |4} A Mathematica notebook for all of the calculations is 
available as an electronic supplement [26J. 

In principle, the worldtube could contain any classical radiating system, so it should 
be kept in mind that the only parts of Section [2] that are truly specific to the Lorentz- 
Dirac example are the values quoted for the coefficients in the two series expansions of 
<9 s P C a P (s, r) (equations (14) and (23)). 

The unit system used in the paper is the SI and the Minkowski metric is taken to 
be f] = diag(— 1, 1, 1, 1). Thus, any equation here can be converted to the geometrized 
Gaussian CGS units used by Hawking and Ellis [27], Misner, Thorne and Wheeler [28] 
and Wald [29], for example, by setting c = G = 1 and 4neo = 1 (hence also fio = An). 



2. The tube method 



We start by supposing the particle has a physically sensible (but unknown) classical 
structure, for example, as a soliton in some nonlinear field theory that reduces to 
Maxwell electrodynamics in its weak- field limit, or perhaps as some classical electron 
model [16]. Here, physically sensible is taken to mean there exists a conserved stress- 
energy tensor T for this underlying theory, and that the particle has a finite size. Thus, 
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we suppose 

V-T = 0, (2) 




rp I matter fo r T ^ ?"0 j /o\ 

' otherwise , 

where T e . m . is a purely electromagnetic (vacuum) stress-energy tensor and the radius 
r defines, for some predetermined level of accuracy, the extent of the matter fields or 
size of the particle structure. A precise value for this structure radius is not needed. It 
will suffice to assume that such an r exists, so that the particle can be considered to 
lie within a worldtube of radius r . The worldline C on which this tube is centred will 
serve as a reference worldline for the particle, about which the asymptotic field of the 
particle shall be prescribed, usually in terms of a multipole expansion [TUl H51 EOl EH] ■ 

The simplest such prescription is that the charged particle have an asymptotic 
field that exactly coincides with that of an electric monopole, as given by the retarded 
Lienard-Wiechert potential for a point-source with worldline C. Later, we shall see 
that this case leads directly to the Lorentz-Dirac equation 0. More generally, the 
asymptotic field of the (non-zero size) charged particle would include higher order 
multipoles, and the corresponding self- force would differ from the Lorentz-Dirac result. 

The worldtube for which Dirac's balance relation is to be derived will be chosen to 
have a tube radius r that is much larger than the structure radius r , so that it lies in 
the asymptotic field of the particle, where the prescribed multipole expansion can be 
used. On the other hand, we still take r to be small compared with the scale set by 
the radius of curvature of C, that is \rnx\ <C 1, where ki (defined by (38)) is the first 
Frenet-Serret curvature of C. The situation described is shown in Figure [TJ 

Let C be a smooth timelike worldline with equation x = z(s) where x a , a = 0, ... ,3, 
are Minkowski coordinates with metric rj = diag(— 1, 1, 1, 1), and the parameter s is 
proper distance (c x proper time) along C, so that z ■ z = —1. In Section [3j spherical 
polar tube coordinates {s, r, ip} are constructed using the Frenet frame of C, such 
that the worldline parameter s coincides on C with the spacetime coordinate s, and 
C has the equation r = 0. The various integration surfaces that are needed can then 
be specified in terms of coordinate surfaces. The worldtube S tube (s,r), the tube-cap 
^cap^r), their intersection 2-sphere Q(s,r), and the 4-volume V(s,r) are defined by 

£ tube (s,f) = {s e [0, s], r = f}, (4) 

£ cap (s,f) = {s = s, r 6 [0, f]} , (5) 

Q(s, f) = {s = s, r = f} , (6) 

V(s,f) = {se[0,s},re[0,f}}, (7) 

where in each case the angular coordinates assume the range ($,<p) G [0, it] x [0, 2ir). 
Figure [T] shows these various structures for the time-symmetric tube coordinate system 



defined in Section |3.1.2| For this coordinate system the worldtube is of Dirac-type, 
meaning that the tube caps are orthogonal to C, as used in pQ. The calculations can 
also, for example, be based on retarded (advanced) tube coordinates, in which case the 
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Figure 1. The worldtube £tube(s,f) used in the self-force calculations. The tube is 
of length s and has radius r > ro where ro is the structure radius of the particle. The 
momentum P C ap(s,r) includes unknown contributions from the stress-energy tensor 
associated with the particle structure. Nevertheless, expressions for the derivatives 
d s P tu hc(s, r) and <9 r P cap (s, r) can still be evaluated in the asymptotic region r ^> r , 
and the self-force on the radiating particle thereby calculated using conservation of 
energy- momentum over the 4- volume V(s,r). 



tube caps are parts of forward (backward) null cones, as described in Section 3.1.3 The 
results quoted in the present section are for the Dirac-type worldtube of Figure [T] 

Let f C ap( s ) r ) be the momentum flux through E cap (s,r), and let Ptube(s,r) be the 
(outward) momentum flux across the worldtube E tube (s,r). In Section 3.2.1 it is shown 
that 

1 



P 



p V' 1 '; '"J 
Ptnbe(s,r) 



Jn(s,r') 



T ■ z/ cap dQ dr' 



T ■ z/ tubc dQ ds' 



(8) 
(9) 



s'=0 JCl(s',r) 



where dfl is the element of area for the 2-sphere Q(r,s), and the vectors z/ cap and 
Vtube (given by (62) and (65)) are normals for the integration 3-surfaces. Observe that 
expressions (|8j) and ^ depend on r and s respectively, only through the upper limits 
of the outer integrals. Differentiating (J8]) with respect to r and ^ with respect to s 
therefore gives 

1 



<9P C ap(s,r) 
dr 



T ■ z/ cap dfi 



(10) 
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ds 



T ■ z/t u bc dfi . 



Following Dirac, equation (11) becomes a relation involving P ca . p (s,r) by considering 
energy-momentum conservation for the fields within the 4- volume V(s,r). Dirac 
assumes these fields are purely electromagnetic, so that T = T e . m . all the way up to 
the worldline C, but in our case only the asymptotic field of the particle is known to us. 
The particle structure is unknown, as is the stress-energy tensor T matter associated with 
this structure. Nevertheless, Gauss' theorem still applies and in Section 3J2 it is shown 
that by integrating ^ over the 4-volume V(s,r) one has 

Pr»Js,r) = P C ap(0,r)-P tube (s,r). (12) 



cap 



Differentiating (|12|) with respect to s and making use of (11) then gives 

1 



<9P cap (s,r) 
ds 



T ■ v tuhc dtt 



(13) 



Si(s,r) 



We now use the fact that the worldtube has radius r 3> r . The 2-sphere Q(s, r) over 
which the integral (13) is evaluated therefore lies in the asymptotic field of the particle, 
where the stress-energy tensor is of known electromagnetic form. After replacing T 
in (13) by T e . m ., one then expands the integrand T e m • v t ube, as a power series in the 
tube radius r, which although large with respect to the structure radius r , is still 
assumed to be small with respect to the scale set by the worldline curvature (\rtz\\ <C 1). 
The angular integrations then involve trigonometric polynomials in if} and may be 
evaluated exactly. In Section [4] these calculations are described in further detail for the 
Lorentz-Dirac example, in which taken to be the stress-energy tensor of an 

electromagnetic field that is the sum of a background field F and an electric monopole 
(retarded Lienard-Wiechert) field for the particle. The result is Dirac's balance relation 
(equation (18) in PQ), 



9P cav (s,r) 
ds 



„ CO 

/ T e . m . ■ z/t u be dQ = ^ a k (s)r l 

Jn(s,r) ,_ i 



>n(s,r) k __ 

where the first few coefficients in the series are found to be 

KiJVi, 



a-i(s) 
a {s) 
aUs) 



87re c 
qF ■ N + 



[k l N l + k 1 k 2 N 2 ] 



Q 



67reoc 



48ne c 

— 4 + 2n 2 k 1 ) N 2 — Ak 1 k, 2 k 3 N 3 



(14) 

(15) 
(16) 

(17) 



Here and K a are respectively the Frenet-Serret frame vectors and curvatures of C 
(see Section 3.1.1 for definitions). For k > 2, the coefficients a^(s) involve successively 
higher derivatives of F evaluated on C, and become rapidly more complicated. 
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A similar calculation is used to evaluate (10). Again, on noting that Q(s,r) lies 

expand the integrand as a 



in the asymptotic field of the particle we replace T by T e 
power series in r and then do the angular integrations. The result is 

<9fcap (■§,?") 

9r c j n ( Sjf .) 



If °° 

- / T e . m . • z/ cap dtt = ^ b k {s) r k , 



where the first few coefficients in the series are 

N , 



b-i(s) 
b (s) 



87re c 
0, 

q 2 



(18) 

(19) 
(20) 
(21) 



487re c 

Assuming convergence of the series, we now integrate ( [l~8| as a partial differential 
equation in independent variables {s,r} to get, 



-Pca P (s,r) 



p(s) 



~2{s r 



" l + Y J bk-i{s)r k /k. 



(22) 



k=l 



where the quantity p(s) (independent of r) appears as the integration "constant". It 
turns out that if the tube caps S cap (s,r) are orthogonal to the worldline C, as they are 
for the Dirac-type tube being used here, then this integration constant can be directly 
identified with the momentum of the particle (the situation for a worldtube with null- 



cone caps is covered in Section 2.2). Differentiating (|22j) with respect to s finally gives 

dP ca ,p(s,r) ^ ^ £ , i. /„\„Jfc/7. ^23) 



ds 



p(s)-6_ 2 (s)r 1 + J2bk-i(s)r k /k. 



k=l 



We have thus found a second expression for the quantity <9 s P cap (s, r), in addition to 



Dirac's balance relation (14). Subtracting (23) from (TT4J) then gives 







a_i + 6_ 



,-1 



(o 



OO 

p) + f - ^-lA 



(24) 



fe=i 



Since r is arbitrary, vanishing of this series entails that all coefficients in the series are 
zero. The r°-term gives an equation to be satisfied by p(s), whereas the terms involving 
r vanish identicall}|§J Using the Mathematica code one can check this is so to at 
least 0(r 5 ). We remark that to explicitly verify that the coefficient of r k is identically 
zero requires, for k > 2, that on C the background field F must satisfy the source-free 
Maxwell equations and also derivatives of these equations (up to a derivative order that 



increases with k). The vanishing of the r°-term in (24) gives the momentum balance 
equation, p = a . For the Lorentz-Dirac example, the coefficient a is given by (16) and 
this equation becomes 



p{s) 



qF-N + 



{k 1 N 1 + k x k 2 N 2 



(25) 



67reoc 

For k 7^ 0, the coefficients etfc(s) are therefore perfect differentials. Having noticed this before 



obtaining (24), the author wrote a Mathematica package for symbolic integration of Frenet frame 



expressions such as ( 17 1 . This package is available as part of 
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This is an exact statement of momentum balance for the case in which the asymptotic 
field of the charged particle is equal to the retarded field of an electric monopole. 



Note that p(s) is still unspecified. To obtain an equation of motion from (25), one 
must specify how the particle momentum is to be defined in terms of the geometry of 
the worldline and the internal degrees of freedom of the radiating system or particle. 
The Lorentz-Dirac equation is obtained by assuming (as does Dirac) that the particle 
has no internal degrees of freedom and that p(s) takes the simplest possible form, 

p(s) = cmN (s). (26) 

Then (25) becomes the Lorentz-Dirac equation, the equivalence of (25) and Q being 



easily established using the Frenet-Serret equations (38)-(41). 

There are several logical gaps in the above derivation of the Lorentz-Dirac equation, 
but these gaps are not related to evaluating the self-force. If one knows the asymptotic 
electromagnetic field of the extended particle (even if this field is defined so as to depend 
on internal degrees of freedom of the particle) then the tube method can be used to derive 



an exact statement of momentum balance analogous to (25), that defines the self-force. 
The gaps only become evident when attempting to derive an equation of motion from 
this balance equation. For example, definitions of the momentum and center-of-mass 
worldline of an extended particle have been given by Dixon in [32] [S]. But our p(s) 



is defined as the the r- independent term (the integration constant) in the series (22) 
for -P cap (s, r). Is our p(s) consistent with Dixon's definition of momentum? Does the 
worldline C have to be chosen to coincide with the center-of-mass worldline? We leave 
these questions for future investigation. For more on the problem of deriving equations 
motion for radiating extended charged particles, see [35J. 

2.1. Classical mass renormalization 

The classical mass renormalization procedure can be understood as a short-cut that 



amounts to assuming an O(r) approximation to our expression (22). To see this, we 



reproduce here the argument used in pQ to get from Dirac's balance relation (14) to the 



result (|25j). 

The traditional interpretation of the quantity P ca , v (s,r) has been that in the limit 
r — > this represents a "bare momentum" of the particle. Assume now that as r — > 
the bare momentum takes an almost familiar form, being the product of the 4-velocity 
of the particle with a bare mass function m g (r). Thus, (compare equation (20) in [Tj) 
let 

P cap (s, r) = cm B (r)N (s) + O(r) . (27) 

One now chooses the bare mass to be negative and diverging as r — > in exactly the 
manner required so as to cancel the problematic r~ l term appearing on the right-hand 



side of (14). Moreover, the difference between these divergent terms must account for 
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the physically observable mass of the particle. Thus, (compare equation (21) in pQ) one 
postulates that 

.2 1 



+ cm , 
87reoC r 

where m is now referred to as the "renormalized" mass of the particle. 
((281) into ((271) gives 



^ca P (s,r) = cmN Q - 



Q 



8ne c r 



+ 0(r) 



(28) 
Substituting 

(29) 



which is recognised as being the O(r) approximation to (22), for a particle momentum 
given by (26). So an O(r) approximation to the Lorentz-Dirac equation is obtained 

The remaining r-dependent terms are then 



0, in which case (28) becomes 



by following the steps previously given 
eliminated from the equation of motion by taking r 
an infinite mass renormalisation. 

The renormalization procedure evidently depends on being able to antidifferentiate 
with respect to s any singular terms that appear on the right hand side of (13) so that 
these can be absorbed into a redefined momentum on the left hand side of (13). In 
1944, Bhabha and Harish-Chandra [I] showed that for point particles possessing any 
multipole moments whatsoever, this can always be done. Using our present notation, 
what they found is that (compare equation (19) in |4]), 



d_ 
Or 



T ■ v tuhc dVt + 



U(s,r) 



d_ 
ds 



T ■ z/ cap dVL 







(30) 



which is the relation (d r d s — d s d r )P cap (s, 



expressed in terms of (10) and (13). 



They go on to observe that if T can be expanded as a power series in r, then (30) shows 
that, with the exception of the r°-term, all coefficients in the series for (13) (which 
they call "the inflow") are perfect differentials. In a later paper [5], the authors refer 
to this result as their inflow theorem. In particular, in [lj they say "the singular terms 
being perfect differentials can always be compensated by the addition of suitable terms 
to P cap (s, r )"[]}] These two papers are both concerned with point particle theory (the 
phrase "point particles" even appears in both titles) and there is no indication that the 
authors were aware of the full significance of their inflow theorem: If all of the perfect 
differentials, not just the singular terms, are "compensated" then all dependence on r 
vanishes, exactly as it did in equation (24). There is then no reason to take r — ► 
and the resulting balance equation (the counterpart of (25)) is exact for finite radius 
worldtubes. Reference jl] has received little attention (an exception is [10] page 361, 
but again, in the context of renormalization for point particles) and the fact that the 
inflow theorem underlies a finite radius tube method seems to have been overlooked. 



2.2. Identification of the particle momentum 



It was remarked earlier that being able to identify the integration constant in (22) as 
the particle momentum depended on having chosen to use a worldtube with caps that 



The original text refers to a quantity A(t) that corresponds to our -P cap (s,r) in (13 1. 
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are orthogonal to the worldline. To see how a result that is independent of tube cap 
geometry can come about, one can do the calculations using a variety of worldtubes, with 
differing tube cap geometries, then demand that the result not depend on any particular 
choice of worldtube. Thus, in Section [3] we define a 1-parameter family of worldtubes, 
with parameter e G [—1, 1]. The tube caps vary with e, interpolating between forward 
null cones for e = 1 and backward null cones for e — — 1. The Dirac-type worldtube, 
with caps orthogonal to C, corresponds to e — 0. 



In this more general setting, the integration constant in (22) depends on e, and 
thus p(s) in (22) is to be replaced by p £ (s). In place of (25) one then finds that |26j, 

p £ (s) = qF-N + (-e Kl 2 N + (1 - e) {k.N, + Ki k 2 N 2 )) . (31) 

An equation that is independent of e is obtained by integrating the ^-dependent terms 
and absorbing them into the definition of the integration constant. Thus, by setting 

,2 



Pe[S) 



p(s) 



q 



6ti€ c 



£K 1 N 1 



(32) 



one recovers the relation (25). 



2.3. A sanity check — the weak energy condition 



In Section |2.1| we saw that if the classical mass renormalization procedure is fully 
implemented using Bhabha and Harish- Chandra's inflow theorem, then one recovers 
(in a somewhat awkward manner) the finite radius tube method of Section [2j Taking 
r — > is then quite pointless and only serves to make the theory mathematically ill- 
defined and physically nonsensical. In particular, the bizarre notion of negative bare 



mass (equations (28)-(29)) arises because the worldtube is too small to fit the particle 
within it. Using the weak energy condition, this can be turned around to determine 
a minimum possible size for the extended particle or system by finding the smallest 
worldtube that could enclose a system with the same prescribed multipole properties. 

If at every point the stress-energy tensor obeys the inequality T af3 u a u 13 > for 
any timelike vector u, then T is said to satisfy the weak energy condition [27]. This 
is equivalent to the energy density being non-negative for any observer. Assuming the 
weak energy condition for T, one has from ([8]) that 



N -P cap (s,r) 



r'=0 Jn(s,r') 



T af3 N a N^dndr' 



> 0. 



where we have used a Dirac-type worldtube, for which (62) gives v^ 



'cap 



(33) 



-iV n 



Substituting expression (22) for P C ap, which we saw in Section 2.2 is specifically for 



a Dirac-type worldtube, one has 

- N -p + N •6_ 2 (s)r- 1 



O(r) > 0. 



(34) 



If the asymptotic field of the particle is that of an electric monopole, then b-2(s) is given 



by (19). In this case, after dropping O(r) terms and rearranging, (34) can be written as 



r > 



cm 



iireo mc 2 (—N ■ p) 



(35) 
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It is interesting to see what this result says about classical models for the electron. In 
the Lorentz-Dirac model the momentum is assumed to be tangent to the worldline of 
the particle and thus, cm/(—N ■ p) = 1. Inequality (35) then becomes r > r e /2, where 
r e is the so-called classical electron radius |36j, 

1 e 2 

r e = = « 2.8 x l(r 15 m, (36) 

47T€o m e c l 

where q = —e is the electron charge and m = m e is the electron mass. However, 
it is known experimentally [37] that the minimum tube radius for an electron is at 
least 10 3 times smaller than r e /2. Not surprisingly, it follows that the Lorentz-Dirac 
classical electron model would be completely inadequate for describing the physics in 
such experiments. On the other hand, our result (35) can not be used to rule out all 
classical models of a spinning electron. In a classical spin model the worldline is typically 
a spacetime helix. The particle momentum is not tangent to the worldline, but instead, 
is directed along the helix axis. The factor cm/ (— iV -p) appearing in (35) can therefore 
be arbitrarily small if the circular motion of the charge centre is close to the speed of 
lighij^j] Such models arise naturally in 2nd and higher order Lagrangian mechanics (for 
example, see [38| 139]). 



3. Geometry of the worldtube 

At least two obvious choices present themselves for the worldtube. One is a worldtube 
with caps orthogonal to C, as used by Dirac |1], the other is a worldtube with caps 
that are parts of forward null-cones based on C, sometimes called a Bhabha tube and 
used, for example, in jl], [12] and [TT] . In either case, E tu b e (s,r) and E cap (s,r) can be 
defined by equations (HJ) and ^ as coordinate surfaces in a system of spherical polar 
tube coordinates {s, r, tp} centred on the worldline C, with the angles <p} defined 
with respect to the Frenet frame of C. 



3.1. Spherical polar tube coordinates 

3.1.1. The Frenet frame of C Recall that the world-line C has the equation x = z(s), 
where x a are Minkowski coordinates and s is proper distance along C. Let the Frenet 
frame vectors of C be 

N^s) = N/(s)e p , (37) 

where e^ are Minkowski basis vectors, e a ■ e^ = rj a/3 . The defining relations for the 
Frenet frame are that iV = z = dz/ds is the future-directed unit tangent to C; that the 
frame N i is orthonormal and has the same orientation as the Minkowski frame e^; and 
that the frame satisfies the Frenet-Serret equations 

N = KiJVj , (38) 

Calculating in the momentum rest frame, one has p = (cm, 0) and N = 7(1, v/c) where 
7 = (1 — (v/c) 2 ) -1 / 2 and v is the 3-velocity of the circular motion. Then, cm/(—N ■ p) — 7 _1 . 
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Ni = KiiVo + K 2 N 2 , (39) 
N 2 = —k 2 N 1 + k 3 N 3 , (40) 
N 3 = -k 3 N 2 . (41) 

The scalar fields K a (s) are the Frenet-Serret curvatures of C. Differentiation with respect 
to proper distance s is indicated by an overdot. The orthogonality conditions for the 
Frenet frame are N i • = rj^, where T] = diag(— 1, 1, 1, 1), or in terms of components, 
N i a N ; j l3 r] al3 = r]^. The matrix inverse of is therefore given by its Minkowski 
transpose, N k a = r\ kj 'Nj l3 r] af3 . Finally, one has the following relation between the 
Minkowski basis covectors and the Frenet coframe, 

e a = WNi" , (42) 

where e a = T] a ^e^, and N l = rj^Nj . 

3.1.2. Time-symmetric tube coordinates Time-symmetric tube coordinates {s, r, ip}, 
centred on the world-line C are defined by the coordinate transformation 

x = z(s) + r n(s, (p) . (43) 

Here x = x a e a and z(s) = z a (s)e a are Minkowski position vectors, and n is a unit 
vector orthogonal to C, parametrized in terms of polar angles as 

n = N^s) sin $ cosy? + N 2 (s) sin $ sin ip + N 3 (s) cos $ . (44) 

The surfaces s = const, are 3-planes orthogonal to C at x = z(s). Within these 3- 
planes the coordinates {r, tp} are standard spherical polar coordinates, with the polar 
angles measured with respect to the directions defined by the Frenet frame at x = z(s). 
We call these coordinates "time-symmetric" rather than "instantaneous rest-frame" or 
similar because the Frenet frame is Fermi-propagated (non-rotating) along C only if 



= 0. Note also, that implicit in M3J) is the identification (s, r, cp) = A^ (s) 



that is, the frame at x — x(s,r, ^, <p) is obtained by parallel transport of N i (s) from 
the point x = z(s). Thus, the Frenet frame, originally defined only on C, is henceforth 
considered to be a frame field on Minkowski spacetim^J 

3.1.3. Retarded (advanced) tube coordinates Retarded (advanced) tube coordinates are 
obtained by taking the s = const, surfaces to be the forward (backward) null-cones based 
on C. The coordinates {s, r, <p} are defined by the transformation 

x = z(s) + r(n(s, 0, <p) + eN (s)) , (45) 

where e = 1 gives the retarded coordinate system and e = — 1 gives the advanced 



system, as shown in Figure 2. Transformation (45) can also be interpreted as defining 



a continuous 1-parameter family of tube coordinate systems. Since e = corresponds 



to the time-symmetric tube coordinates (43), we may deal with all interesting cases at 



once by carrying e 6 [-1,1] as a parameter through our calculations. 

+ This frame field is only needed in the vicinity of C, where the tube coordinate metric g is non- 



singular. For the coordinate transformation (45), which includes (43) for e = 0, one has yj— det(g) 



' sini? (1 + r^ 1 (l — e ) sini9cos^j) near C, so g is non-singular if \rKi\ < (1 — e )~ 
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Figure 2. The tube caps are s = const, surfaces in one of the tube coordinate systems 
defined by the transformation x = z + r(n + sN ). Shown here are the null caps for 
e = ±1, and the Dirac-type cap for e = (the r — const, worldtubes in this figure 
coincide for different e only because the figure has been drawn for a source in linear 
motion) . 



3.2. Gauss' formula for the worldtube 



We require Gauss's formula specialised to integrating over the tube volume ([7]). Writing 



the boundary of V(s,r) as dV(s,r) 
for any C 1 vector field X, 



^ cap (s,r) - E cap (0,r) + E tube (s,r) one has that 



V-XdV(s,r) 



V(s,r) 



X -dZ 



cap 



s,r) - 



E ca p(s,r) 



>(0,r) 



X-cffi cap (0,r) 



- / X ■ dE tuhe (s,r) . 

'E tubo (s,r) 



(46) 



In Sections 3.2.2-3.2.4 it is shown that these surface integrals can be evaluated as 



X -dZ r 



X ■ i/ cap dfi dr' , 



(47) 
(48) 



J cap( s ?' r ) — 
Jcap (s,r) Jr'=0 Jn(s,r') 

X ■ dE tube (s, r) = / / X ■ z/ tube dfi ds' . 

'S tube (s,r) Js'=0 Jn(s',r) 

Here dfl is the volume element of the 2-sphere Q(s,r). It follows that the orientations 
of Et u be(s,r) and the final (s > 0) tube cap E cap (s,r), are given by the directions of the 
normals z/ tube and z/ cap respectively, and for (46) to hold these are to be chosen so that 



X ■ v is positive if X is a vector that points out of V(s, r) (see for example Section 2.8 
in [27], or Appendix B.2 in [29]). The initial (s = 0) tube cap is then reverse-oriented 
with respect to ^ cap since X ■ z/ cap is negative for outward pointing vectors on E cap (0, r). 



This is in accordance with the sign of the corresponding term in (46). The normals are 



given by (62) and (65). 



3.2.1. Energy-momentum conservation Let T = T al3 e a <S> be the conserved stress- 
energy tensor (J3J) . The four vector fields (for a = 0, . . . , 3) defined by X^ = e a ■ T = 
T^e^ are then divergence free, V ■ X^ = 0. Applying Gauss's formula (46) to each of 
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these vector fields gives the four conservation equations, 
= P c %(s,r)-P" p (0,r)+P« he (s,r), 



where 



P c t P (s,r) = - [ X^-dS cap ( S ,r) 



P 



C JS cap (s,r) 

,nj,eis,r) = - J X {a) • dE tube (s, r) . 



e « P cap and Ptube 



(49) 



(50) 



(51) 



e a J tube- 



The momentum flux 4-vectors are defined by P cap 
Equations (|8])-([9]) then follow from (47)-(48), and the energy-momentum conservation 



law (12) follows from (49). 



3.2.2. The directed hypersurface element dE Suppose we are given a hypersurface E 
defined in terms of spherical polar tube coordinates. That is, E is parametrized by 
coordinates y A , A = 1,2,3 and has equations of the form (s,r, $,ip) = (f a (y A )). We 
shall derive here the Frenet frame expression for the directed surface element of E. 

In Minkowski coordinates the directed hypersurface element is dE = e M dE M , where 
the 3-forms dE^ are given by 

1 



, 7 dx a A da;^ A da; 7 . (52) 

Here e p is the fully antisymmetric Levi-Civita tensor (£ 012 3 = 1 in Minkowski 
coordinates). For E parametrized by coordinates y A , the 1-forms da; Q in (52) can be 
evaluated as 

° (e a ■ x) dy A = e a .^dy A . (53) 



dx a 



dy J 



dy J 



The vectors dx/dy A that appear here are to be calculated by differentiating (45), so are 
naturally obtained as linear combinations of the Frenet frame vectors. Let these linear 
combinations be 

(54) 



Using (42), equation (53) now becomes 

dx a = e a -N„X\dy A = N,: Q X j A dy J 



(55) 



Again using (42), one has dE = iVW^dE^. Substituting (55) into (52) then gives 

dE = ^N i e mf3l N/N a a N b l3 N ( ;X a A X b B X c G dy A Ady B Ady c , 

and because the Frenet frame is orthonormal and has the same orientation as the 
Minkowski frame, 

Wr^/V^c 7 = e iabc . (56) 
The hypersurface element therefore simplifies to 

dE = N l e iabc X\X b 2 X% dy 1 A dy 2 A dy 3 , (57) 
which is the required Frenet frame expression. 
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3.2.3. Directed surface element for a tube cap The tube cap E cap (s,r) is denned by 
(|5| as part of the s = const, coordinate surface. For e = the caps are orthogonal to 
C, whereas for s = ±1 they are null-cones based on C. As parameters for the tube cap, 



one may take {y A } = {r, $,<p}- The derivatives dx/dy A are found using (45), 

dx/dr =n + eN , (58) 

dx/dd = r (N l cos i? cos ip + N 2 cos d sin <p — N 3 sin i?) , (59) 

dx/dp = r (—N l sin •& sin p + N 2 sin $ cos <£>) . (60) 



From these expressions one can read off the coefficients X 3 A , as defined by (54). The 



directed hypersurface element (57) is found to be 
d£ 



z/ cap dfi A dr 



where dfl = r 2 sin •& dd A dtp, and the normal is 



'cap 



-(Wo + en) 



(61) 



(62) 



where the sign of u cap has been chosen as described in Section (3.2). 



3.2.4- Directed surface element for the worldtube The worldtube S tubc (s,r) is defined 
by p]) as part of the r = const, coordinate surface. As parameters for the worldtube, 
one may take {y A } = {s, <p}. The derivatives dx/dy A that are needed to evaluate the 



surface element are (59)-(60) together with 



dx/ds = N + r (h + en l N l 



(63) 



where 



n = (k-^Nq + K 2 iV 2 ) sin •& cos <p + (— k 2 N x + k 3 N 3 ) sin -dsimp — k 3 N 2 cos . 
For the coordinate ordering specified by {y A } = {s, ■$,<£>}, one finds that (57) gives 



dS = — z/ tu bc dfi A ds . 
where the normal 

^tube = n + tk x sin i? cos (n + eA^ ] 



(64) 
(65) 



has its sign chosen as described in Section (3.2). 



4. The stress-energy tensor 

For r > r , equation ([3]) defines T to be the stress-energy tensor of an electromagnetic 
field. We take this field to be the sum of a background field F and a prescribed retarded 
field .F ret for the particle or radiating source that we are interested in, 

F tot = F + F ret . (66) 

The background field F will include any applied field, any fields generated by other 
moving charged particles, and also any radiation from the source that has been reflected 
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back to the source positiorj^jby boundaries or material media. Thus, F will appear in the 
calculations as an unspecified solution of the Maxwell equations, known to be source-free 
throughout the the spacetime region occupied by the particle or system of interest. 

The stress-energy tensor corresponding to F tot is T e . m . = T£j m N i <g> iV- , where the 
frame components are 

TiL. = ^ (F&F* Vat ~ \ rpF&Ft&Wcd) • (67) 

For the Lorentz-Dirac charged particle, the field of the radiating source is taken to be 
the retarded Lienard-Wiechert field. 



4-1. The Lienard-Wiechert field 

The retarded Lienard-Wiechert potential for the electromagnetic field of a point charge 
q with worldline C is 

M*) = J 3 — %# , (68) 
4ire c R(x) 

where the vector and scalar fields z Tet (x) and R(x) are defined as follows. Let the 
retarded proper distance s ret (x) be the scalar field that gives the worldline parameter 
value of the source-point on C for a retarded field at x. Thus, s Tet (x) is defined so that 
x lies on the forward null-cone based at z(s mt (x)), and 

C0&) = x- z(s iet (x)) , (69) 

is therefore a future directed null vector field. Then z rct (x) = z(s ret (x)) is defined as the 
unit tangent to C parallel propagated from the retarded source point z Tet (x) = z(s Tet (x)). 
Similarly, one defines z Tet (x) = z(s iet (x)). The scalar field R is the radial distance of 
x from the source, as measured in the instantaneous rest-frame of the source at the 
source-point z iet , 

R = -iret • C • (70) 

The scalar and 3- vector potentials are defined by A a = (<&/c, A). For a charge at rest, 
i r g t = (1,0,0,0), and (68) defines $ to be the Coulomb potential. 
The electromagnetic field tensor corresponding to (68) is 

•^retajS = Ap^ — A^p. (71) 

The coordinate derivatives of the potential can be calculated using the relation 

Srct(x) = -QR' 1 = -kp , (72) 

that follows from differentiating the equation Q^CVav = with respect to x 13 . Here we 
have introduced the null vector k = (R~ x , satisfying the normalization k ■ z rct = —1. 

* This last contribution to the background field accounts for changes in spontaneous emission rates 
and effects thereof, such as bandgaps in photonic crystals. 
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ij 



The Frenet frame components of the field tensor are then given by F ret 
One finds that 
Q 



r ret 



((& y - z r { t k*) (1 + aR) R- 2 + (4 t k j - zi t k l ) R' 1 ) 



where a is the acceleration dependent scalar field 

O = ^ret ' • 

The frame components of the derivative fields i ret 



(73) 



(74) 



^ret 



iet Nj can be 



read from the series expressions (76) and (77) of the following section. Likewise, £ is 



given by the series ( 75 ) and the scalars R and a then follow from ( 70 ) and ( 74 ) 



4-2. Tube coordinate expressions for retarded fields 

Recall that the retarded proper distance s iet (x), is such that x lies on the forward null- 
cone based at z(s rct (x)). The retarded tube coordinate system is constructed so that 
Sret(^) — s - In an Y other tube coordinate system one must solve for s ret (x) using the 



condition that ((x), given by (69), is a null vector field 



For points sufficiently close to C, the separation between s and s Ict (x) will be small. 



Thus, making use of (45) and writing s ret (x) = s + A, one expands (69) as 
( = r{n + eN )- (z A + z A 2 /2! + z A 3 /3!) + 0(A 4 ) 



A 2 
2! 



r (n + eN ) - N A - N, k x =p - (JV ^ 2 + ^ « x + iV 2 % « 2 ) ^ + 0( A 4 ) . (75) 



Let A = YHt=i c k rk i f° r coefficients Ck = Ck(s,$,<p) to be determined. Substituting the 



series for A into (75) gives £ as a series in r. Then, by forming the series for the null 



expression (■( and equating coefficients to zero, one obtains a sequence of equations that 
can be solved in turn for the unknown coefficients The first coefficient is found to be 
c\ = e±l. Since the proper distance parameter has been taken to increase to the future, 
the two solutions for c\ give s iet (x) = s— r(l— e)+0(r 2 ) and s a d v (a ; ) = s+r(l+e)+0(r 2 ). 
That is, the retarded solution corresponds to c\ = e — 1. For brevity, let fi = e — 1. The 
first few coefficients in the series for A are then 

ci = /i, C2 = - /i 2 K 1 sint? cos 99 , 

C3 = — (— (4 + 3/i) K]^ 2 + 3 (4 + /i) (k 1 cos ip sin 1?) + 4 sin ■§ (k 1 k 2 sin <£> + % cos </?)) . 

Observe that in retarded tube coordinates (for which \i — e — 1 = 0), one has = 
and thus expected. 

Using the above solution for s iet (x) = s + A, the fields z Tet (x) and £ re t(x), needed 



in expression (73) for the frame components of the retarded electromagnetic field, can 



be obtained by expansion in r. One finds that 

U 2 

z mt (x) = N + ^N x K x r + y {N k 2 + N x (z^ 2 sintfcos^ + >%) + N 2 k x k 2 ) r 2 + 0(r 3 



(76) 



^ret 



N 1 k 1 + n (N k 2 + N 1 k 1 + N 2 k x k 2 ) 
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+ — ^A* r K 1 (fi^ 2 sin i? cos ip + 3/%) + iVj (fi^ 3 — k 1 k 2 2 + sin i? cos </? + k^) 



+ iV 2 (/% 2 /«2 sin $ cos if + 2k 1 K 2 + + ^3 K i K 2 K 3j r2 + 0( r 

The resulting expressions for F^ t can be found in [26J. 



(77) 



4-2.1. Electromagnetic field in retarded tube coordinates The field F ret is particularly 
simple in retarded tube coordinates because in this case 

i ret = N , z ret = n 1 N 1 , k = n + N , R = r , a = k x sin 1? cos . (78) 

The frame components of an electromagnetic field tensor are identified with the electric 
and magnetic fields in that frame by 

/ Et/c E 2 /c E 3 /c \ 

-Ei/c B 3 -B 2 

-E 2 /c -B 3 B 1 

\ -E 3 /c B 2 -B l J 

for example, F 01 = E\jc. The electric and magnetic fields for F ret with respect to the 



(F ij ) 



(79) 



Frenet frame defined by the retarded tube coordinate system are easily found from ( 73 ) 



using (78) 



E 
B 



Q 



4:7160 

Q 



y + (n sin •& cos <p — N x ) — J , 



iV 9 cos d + Nn sin ■& sin 09 
47re c V d V / 



(80) 
(81) 



The magnetic field in the frame associated with time-symmetric tube coordinates is 



non-singular. The singular field (81 ) can be traced to mixing of E and B by a boost, of 



approximate magnitude r/t 1; that relates these two different frame fields. 



4-3. Expansion of the background field near C 

The values of the background field F = F a ^e a ® = F^N i ®Nj on the worldtube cross- 
sections fl(r, s) are to be approximated in terms of frame components of the Minkowski 
derivatives of F(x) at the point z(s) on the worldline C. By Taylor expansion, 



i^'(s,r,^) = [F ij ] r=Q + 



dr 



2 r^2 



+ 



r=0 



2! 



«9 2 i™ 
dr 2 



+ ... 



r=0 



The frame field N^s) is independent of r, so the radial derivatives can be calculated as 



dr 



dr dx^ 
l^N^F^N^ 
dx k 



(82) 
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where x k and F %3 k are respectively the frame components of the Minkowski position 
vector N k , and the first Minkowski coordinate derivative of F. Similarly, 

for example, 

d 2 F lj dx k dx l ij d 2 x k {j 
Or 2 = lh~dr~ kl + ~dr^ k ' ( ' 



From (45) one sees that dx k /dr = N k ■ (n + eN ) = n k + e5q, where the components 
of n(s, ip) are given by (44). The second and higher r-derivatives of x k are zero. The 



required Taylor expansion is therefore 



,j -. 2 



F^(s, r, 0, <p) = F%) + r{n k + e 5 k ) F\(s) + - {n k + e 5 k ) (n l + e S l ) F\ (s) + . . . (84) 
where field values on C are denoted, for example, as F^(s) = [F l i(s,r,'&,(p)] r= Q = 

J x=z(s) ' 



N^N^is) [F^(.r 



5. Conclusion 

The paper describes a tube method for calculating the self-force on an extended radiating 
system. The new method is presented for the special case of electromagnetic self-force in 
Minkowski spacetime, but it seems likely that similar techniques could be used to derive 
a tube method applicable to a radiating system of finite size in a curved background 
spacetime. In particular, a finite radius tube method may offer a mathematically sound 
alternative to having to invoke a mass renormalization in the tube-based gravitational 
self-force calculations of [20]. In contrast to the classical mass renormalisation procedure, 
all steps in our tube method are physically well motivated and could be justified to any 
required degree of mathematical rigour. 

For finite size electromagnetic systems in Minkowski spacetime, we find that if 
the field outside the radiating system is known exactly (for example, as a prescribed 
multipole expansion) then our tube method gives an exact expression for the self- 
force. This self-force coincides with that derived by the classical mass renormalisation 
procedure, so in this sense, it has been proved that the renormalisation procedure gives 
correct results in Minkowski spacetime. However, as discussed in Section 2.1[ the reason 



why classical mass renormalisation works is that the point particle limit (tube radius 
r — > 0) becomes superfluous if renormalisation is implemented in full using the Bhabha 
and Harish-Chandra inflow theorem. In fact, a fully implemented renormalisation 
scheme would be exactly equivalent to our finite radius tube method. Thus, taking 
r — > only serves to introduce mathematical inconsistencies and obscure the underlying 
physics. The classical mass renormalisation procedure is therefore conceptually flawed, 
in so far as it is associated with a point particle limit. 

The finite radius tube method follows from Dirac's balance relation together with 



our expression (22) for the momentum flux through the cap of the worldtube. Assuming 



the weak energy condition for the total energy momentum tensor T, expression (22) can 



also be used to determine the smallest possible size of an extended radiating system 



that has an asymptotic field with a given multipole structure. In Section 2.3 this was 
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done for a particle with an asymptotic field equal to that of an electric monopole. We 
found that if such a particle is governed by the Lorentz-Dirac equation, and if its mass 
and charge are set equal to that of the electron, then this particle is at least a thousand 
times larger than the length scale of any possible electron structure, else the weak energy 
condition is violated. On the other hand, the experimental limits on electron structure 
size are consistent with the weak energy condition as applied to helical worldline models 
of a classical spinning electron. 
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